Limit Theorem for Continuous-Time Quantum Walk 

on the Line 

Norio Konn 

Department of Applied Mathematics, 
Yokohama National University, 
79-5 Tokiwadai, Yokohama, 
240-8501, Japan 

(Dated: February 1, 2008) 

Abstract 

Concerning a discrete-time quantum walk x['^'^ with a symmetric distribution on the hue, whose 
evolution is described by the Hadamard transformation, it was proved by the author that the 
following weak limit theorem holds: x['^^ /t dx/iT{l — x'^)\^l — as i ^ oo. The present paper 
shows that a similar type of weak limit theorems is satisfied for a continuous-time quantum walk 
X^'^^ on the line as follows: X^'^^ /t dx/'ir\/l — as f — > oo. These results for quantum walks 
form a striking contrast to the central limit theorem for symmetric discrete- and continuous-time 
classical random walks: Yt/^/i e~^^ ^"^dx/ \/2tt as t ^ oo. The work deals also with issue of the 
relationship between discrete and continuous-time quantum walks. This topic, subject of a long 
debate in the previous literature, is treated within the formalism of matrix representation and the 
limit distributions are exhaustively compared in the two cases. 
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I. INTRODUCTION 



time case 



Quantum walks have recently been introduced and investigated, with the hope that they 
may be useful in constructing new efficient quantum algorithms, (for reviews of quantum 
walks, see Jll 2, Jj)- There are two distinct types of the quantum walk: one is a discrete- 
H H i y, y, H U, Il2|, the other is a continuous-time case H Q Q, Q 
The quantum walk can be considered as a quantum analog of the classical 
random walk. However there are some differences between them. For the discrete-time 
symmetric c/a552ca/ random walk Y° starting from the origin, the central limit theorem shows 
that Y^°/^/t t-'''' 1"^ dx j ^plii as t ^ oo. The same limit theorm holds for the continuous- 



time classical symmetric random walk, (related results can be found in Refs. 20|, |21|). 
The limit density function is the normal distribution and has a bell-shaped curve with one 
peak at the center. On the other hand, concerning a discrete-time quantum walk xf^ 
with a symmetric distribution on the line, whose evolution is described by the Hadamard 
transformation, it was shown by the author j^, |^ that the following weak limit theorem 
holds: Xf^ jt dx/nil - x2)Vl -2x2 as t — > oo. This paper presents that a similar type 
of weak limit theorem is proved for a continuous-time quantum walk XI on the line as 
follows: X'f^/t dx/nVl — as t — s> oo. Both limit density functions for quantum walks 
have two peaks at the two end points of the support. The study of weak limits for discrete- 
time quantum walks is treated also in Ref. |6|, with a simplified proof with respect to earlier 
derivations. 

As a corollary, we have the following result. Let o"^^j(t) (resp. cr'^^^-^it)) be the standard 
deviation of the probability distribution for a discrete-time (resp. continuous-time) classical 
random walk on the line at time t. Similarly, o"^^^(t) (resp. denotes the standard 

deviation for a discrete-time (resp. continuous-time) quantum walk. Then, the central limit 
theorem implies that cr^^^(t), cr^^^(t) x i/t, where f{t) x g{t) indicates that f{t)/g(t) c^{^ 
0) as t ^ oo. In contrast, using our limit theorems, it is shown that (T^^^(t), cr|'^-)(t) x t hold. 
That is, the qunatum walks spread over the line faster than the classical walks in the both 
discrete- and continuous-time cases. 

The remainder of this paper is organized as follows. Sect. 11 gives the definition of the 
walk and our results. In Sect. IE, we prove Proposition 1. Sect. IV is devoted to a proof of 
Theorem 1. Conclusion and discussion are given in Sect. V. 
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II. MODEL AND RESULTS 



Let Z be the set of integers. To define the continuous-time qunatum walk on Z, we 
introduce an oo x cxd adjacency matrix of Z denoted by A as follows: 
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The amplitude wave function of the walk at time t, {'^{t)), is defined by 



\^it)) = umm, 



where 



itA/2 



U{t) = e 



Note that U{t) is a unitary matrix. As an initial state, we take 

|vI/(0)) =^(...,0,0,0,1,0,0,0,...), 

where T indicates the transposed operator. Concerning the details of the definitions for 
the continuous-time case, see li,l2^- Let |\&(/c,t)) be an amplitude wave function at 
location k at time t. The probability that the particle is at location k at time t, P{k,t), is 
given by 



P{k,t) = {^{k,t)\^{k,t)). 

As for results on the walk on a circle, i.e., {0, 1, . . . , — 1}, see j^, 

We will obtain an explicit form of U{t) first. Our approach is a direct computation based 
on the oo X oo matrix A without using its eigenvalues and eigenvectors in order to clarify 
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a relation (stated below) between the continuous-time and discrete-time quantum walks, 
(another approach and the same explicit form can be found in Section III C of Ref. 
Let Jk{x) denote the Bessel function of the first kind of order k. As for the Bessel function, 
see Watson 23] and Chapter 4 in Andrews et al. |i24l] . 

Proposition 1. In our setting, we have 



\ 



U{t) 







-3 
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-1 
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+2 




/•. 
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tJi{t) 


^'J2it) 


v^Ut) 


i^Mt) 


i^'Mt) 


-2 
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Mt) 


iJi{t) 


l^J2{t) 


i^Mt) 


i^Mt) 


-1 






iJi{t) 


Mt) 


lJl{t) 


t'Mt) 


^'Mt) 









t'J2{t) 


iJi{t) 


Mt) 


iMt) 


t'Mt) 


+1 




i^Mt) 


^^Mt) 


l^J2{t) 


iMt) 


Mt) 


iMt) 


+2 






iMiit) 




i^Mt) 


iMt) 


Mt) 



V-- '■) 

That is, the (1,171) component ofU{t) is given by J\i^rn\it) . 

From Proposition 1, looking at a column of U{t), we have immediately 
Corollary 1. The amplitude wave function of our model is given by 

\^{t))=^{...,^'Mt),^'Mt),^Mt),Mt),^Mt),^'Mt),^'Mt),■■■), 
that IS, |^(A;,t)) = i^'^^Miit) for any location G Z and time t >0. 

Moreover, noting that J_fc(t) = {-!)'' Jkit) (Eq. (4.5.4) in Ref. [2^), we have 
Corollary 2. The prbability distribution is 

P{k,t) = Jl^{t) = Jl{t), 
for any location k ^ Z and time t > 0. 

In fact, the following result (see Eq. (4.9.5) in Ref. 

oo 

Jo^(t) + 2^42(t) = l 

k=l 



ensures that Y2T=-oo ^) = 1 t > Remark that the distribution is symmetric 

for any time, i.e., P{k,t) = P{~k,t). 

One of the interesting points of the above mentioned results is as follows. Let C be the 
set of complex numbers. Fix a positive integer r. We suppose that a unitary matrix Ur has 
the following form: 



Ur 







-3 
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-1 





+ 1 
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/•• 
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-2 
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W3 
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Wi 
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W-4 


W-3 
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W^l 


Wo 


Wi 


+2 




W_5 


W_4 


W-3 






Wo 



J 



with Wk E C {k E Z) and Ws = for |s| > r. Then No-Go Lemma (see Ref. jlfll |) 
shows that there is the only non-zero w^ with \w^\ = 1; that is, there exists no non-trivial, 
homogeneous finite-range model governed by the Ur- For example, when r = 1, we have 
= I, Wo = Wi = 0", "|wo| = Ij""^-! = Wi = 0", or = l,Wo = w_i = 0". Thus, 

the model has a trivial probability distribution. However, our homogeneous, but infinite- 
range model has a non-trivial probability distribution given by squared Bessel functions 
(Corollary 2). 

An open problem on quantum walks is to clarify a relation between discrete-time and 
continuous-time quantum walks (see Ref. Q], for example). Another interesting point of the 
results is to shed a light on the problem. To explain the reason, we introduce the following 
matrices as in our previous paper Q : 



where we assume that U = Pj + Qj {j = A, B) is a. 2 x 2 unitary matrix. Here we consider 
two types of the discrete-time case; one is A-type, the other is B-type. The precise definition 
is given in Ref. [3]. Then the unitary matrix of the discrete-time quantum walk on the line 







( a Q] 




^0 




and Pb = 




, Qb = 1 








icO) 




^0 


d I 
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is described as 



id) 



O Pj o o o 

Qj O Pj O O 

O Qj O Pj O 

O O Qj O Pj 

O O O Qj O 



\ 



with 



O 




■J 



for j = A and B. The unitary matrix [/(^d) for the discrete-time case corresponds to U{t) 
for our continuous-time case at time t = 1. More generally, Uf^^-^ corresponds to U{n) for 
n = 0, 1, . . .. Once an explicit formula of U{t) is obatined, the difference between continuous 
and discrete walks becomes clear. As we stated before, U{t) has an infinite-range form. On 
the other hand, f/(d) has a finite-range form. Moreover, we see that U(^d) is not homogeneous. 
It is believed that the difference seems to be derived from the fact that discrete quantum 
walk has a coin but continuous quantum walk does not fl| . However, the situation is not so 
simple, since the discrete-time case also does not necessarily need the coin (see Refs. QQ 
for more detailed discussion). 

We define a continuous-time quantum walk on Z by Xj whose probability distribution 
is defined by P{Xt = k) = P{k,t) for any location A; G Z and time t > 0. Note that it 
follows from Corollary 2 that P{k, t) = Jl{t). Then we obtain a new weak limit theorem for 
a contimuous-time quantum walk on the line: 



Theorem 1. //— 1 < a < 6 < 1, then 
P{a < Xt/t <b) 



IT 



dx. 



x^ 



it 



oo 



Note that 



X 



2m 



1 ttVI — X' 



dx 



Tr/2 



sin^'" (p dip 



for m = 1, 2, . . . , where nil = n{n — 2) 5 ■ 3 ■ 1, if n 

n = even. From Theorem 1 and Eq. (PJ), we have 

1,2 



(2m - 1)!! 
(2m)!! ' 

odd, = n{n — 2) 



(1) 

6 ■ 4 ■ 2, if 



Corollary 3. For m 

E{{Xt/tf"^) 



(2m - l)!!/(2m)!!, (t ^ oo). 
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By this corollary, for the standard deviation of our walk, a'^^-^{t), we see that 

a^^)(t)/t ^ 1/v^ = 0.70710 ... , (t-^oo). 
We consider a discrete-time quantum walk xj^^'^ with a symmetric distribution on the 



line, whose evolution is described by the Hadamard transformation (see Ref. 

1 1 1 



3), 



that is. 



U r 



1 -1 



The walk is often called the Hadamard walk. In contrast with 
the Hadamard walk, the following weak limit theorem holds 



the continuous-time case, for 
BBS : 



Theorem 2. If -1/^2 <a<b< 1/V2, th 



en 



P{a< X^f^ /n<h) / dx, (n-^oo). 

As a corollary, we have 



2 = 0.54119..., fn^oo). 



Comparing with the discrete-time case, the scaling in our continuous-time case is same, but 
the limit density function is slightly different. However, both density functions have some 
similar properties, for example, they have two peaks at the two end points of the support. 



III. PROOF OF PROPOSITION 1 



To begin with, A is rewritten as 



A 







-1 
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+2 






/-.. 
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(2) 
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where 



P 




T 




Q 




This expression corresponds to a unitary matrix for a discrete-time quantum walk in which 
the particle maintains its position during each time step, (see j3] for more details). The 
following algebraic relations are useful for some computations as in the case of a discrete- 
time walk: 



= Q^ = O, PT + TP = QT + TQ = J, PQ + QP = T, 



(3) 



where O is 2 x 2 zero matrix and / is 2 x 2 unit matrix. From now on, for simplicity, Eq. 
((21) is written as 

A=[...,0,0,0,0,Q,T,P,0,0,0,0,0,..]. 

A direct computation gives 

A^ = [ 
A^ = [ 
= [, 

It follows by induction that 



.,0,0,0,0,/,2/,/,0,0,0,0,...], 
.,0,0,0,J,4J,6J,4J,J,0,0,0,...], 
.,0,0, /, 6/, 15/, 20/, 15/, 6/, /, O, O, 



A 



2n 



\ n n 4^^") /i*-^"^ 4^^"^ A^- 



(2n) A2n) A2n) 



,4'"\o,o. 



where A 



(2n) 



A 



for any A; = 0, 1, 

^2n+l _ 



- af )/ and 

2n 

, n. On the other hand, by using A"^"'^^ 



a 



(2n) 



(2n+l) 4(2n+l) ^(2n+l) ^(2n+l) ^(2n+l 



A'_;' 



5 ^0 



A^" X A and Eq. ©, we have 

^VZ'«-|-i; /.(2n+l) ^ -i 

,7^]^ 5 • • • 5 ) '-^) • • -J; 



where 



.(2n+l) 
-("+1) 



a(2")Q, A 



(2n+l) 



-(n-l) 
^(2n+l) 



(2n) 



ar^)T+(ar-ar)Q, 4 



(2n) 



+ a 



(2n). 



T, 



(2n+l) 



(2n) 
'-1 



la; '+4'"))T+(ap-ar^)P,... 



(2n) 



(2n)x 

2 J-' 5 • • • ) 



4(2n+l) _ A2n)rr. , (2n) p A2n+1) _ (2n) p 
— "n "n-l-' ) (n+1) ~ "ri 
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The definition of U{t) gives 



U(t) = e**^/2 ^ y i2i_^n 



n=0 



oo /t\2" oo /^2n+l 



n=0 ^ ' n=0 

Tlierefore we obtain U{t) = B{t) + iC{t), wliere 



with 



B{t) = [..., -B-fc, . . . , -B-i, Bo,Bi,. . Bk, . . .], 
C(t) = [. . . , C-k, . . . , C_i, Co, Ci, . . . , Cfc, . . .], 



and 



C-k = CkI + dkQ, Ck = CkI + dkP, {k = 0,l,...), 



„ \2n+l 

Cfc = 

n=0 

do = 0, 



(2n + l)! tV^-(^-l)7 V^-(^ + l) 
Remark that 

, 2m 



see Eq. (4.9.5) in Ref. |2|) Finally, by using Eq. (0)) and the following relation: 

2n\ f 2n A _ /2n + 1 
k)^\k-l)~\ k 

we have the desired conclusion. 



IV. PROOF OF THEOREM 1 

We begin by stating the following result (see page 214 in Ref. [2^): suppose that 
and c are lengths of sides of a triangle and = + 6^ — 2a6cos^. Then 



Jo(c) = ^ Jk{a)Jk{h)e'^^. 



k = — QO 
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If we set t = a = 6 in Eq. (jS)), then 

oo 
fc=— oo 

From Eq. we see that the characteristic function of a continuous-time quantum walk on 
the hne is given by 

oo 

E{e^^^^) = ^ e^'^^Jm = Ut v/2(l-cosO). (7) 

fc=— oo 

First we consider that t is a positive integer case, that is, ni= t) = 1,2,.... By using Eq. 
((Tj), we have 

as n — s> oo. To know a hmit density funtion, we use the following expression of Jo{x) (see 
Eq. (4.9.11) in Ref. H); 

JoiO = ~ / cos(^ simp) dip = — cos(^ siny^) dep. (8) 
Jo Jo 

Taking x = sin (f, we have 

MO = f cos(ex)— =L= dx = f e'«-— =L= dx. (9) 
By using Eq. (jH)), we get 

2 r'^ 



\M0 - < - [ I cos(esin(^) - 1| 

Jo 



From the bounded convergence theorem, we see that the limit Jo(0 is continuous at ^ = 0, 
since cos(.^ sin y^) — > 1 as ^ ^ 0. Therefore, by the continuity theorem (see page 99 in [2^) 
and Eq. Q, we conclude that if n — > 00, then X„/n converges weakly to a random variable 
whose density function is given by l/vrVl -a;2 for x e (-1, 1). That is, if -1 < a < 6 < 1, 
then 



b 



P{a<Xn/n<b) ^ / — ==dx, 

Ja TTVl - X'^ 

as n ^ 00. Next, to deal with values of t that are not integers, we want to show 
\Pia<Xt/t<b)-P{a<Xit]/[t]<b)\ 0, (t ^ 00), 
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where [x] denotes the integer part of x. To do this, we observe that 



|P(a < Xt/t <b)- P{a < < h)\ 

= \P{at <Xt< ht) - P{a[t\ < X[t] < h[t\)\ 



at<k<bt a[t]<fe<fe[t] 



at<k<bt a[t]<k<b[t] 
yat<k<bt a[t]<k<b[t] 

K^{t) + K^{t). 



< 



a[t\<k<b[t\ 



To estimate K -[ {t) and K2{t), we use the following Meissel's second expansion (see Eq. (5) 
in page 228 of Isaj): 



2cot/3 



cos(g(z/,/5)-7r/4), 



(10) 



where P{v,l3) and Q{v^(3) are defined by Eqs. (1) and (2) respectively in page 228 of 
First we consider Ki{t) case. We see that 

K^it) < max{jf„,](t),4](t),|j2,](t)- J2,](t)|} 
< 2xmax{j2,](t),jJ,](t)}. 



(11) 



Without a loss of generality, we consider only </[^t](^) with < c < 1. From Eq. taking 
z/ = [ct], we have 



2 cot _2p 
7r[ct] ^ 



cos2(g([ct],/3)-7r/4). 



(12) 



where sec /? = t/[ct]{> 1). Then, for large t, we get 

2.^.2 2c 

J\ct\V') ^ r ,1 ^ r. -: 



since cot 13 cj\J\ — as t — > cxo. The last inequality impli 



les 



Jf,,](t)^0, (t^oo). 



(13) 



Combining Eq. (fTT| with Eq. (fT^ implies that Kxif) ^ as t ^ cxd. 
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Next we consider K2(t) case. Noting Eq. (fT^. it is enough to show that 



[at]<k<[bt] 



0, 



(t oo). 



As in the case of ifi()f:), fix < c < 1. Define /3i and P2 by sec Pi = t/[ct] and sec P2 = [t]/[ct] 
respectively. Then we have 



cot Pi ~ cot P2 ' 
cot Pi — cot P2 ^ 

tan Pi ~ tan P2 

tan /5i — tan P2 • 
/?! ~ /?2 ~ arccosc, 



t-[t] 



(1 - c2)3/2 



C 
1 



B B {2-VT^^)c t-[t] 

Di — Do ~ , X , 

where /(t) ~ g{t) means /(t)/(7(t) ^ 1 as t — 00. By using the above estimates, Eq. (O, 
and noting P{\ct],Pi) ^ (t ^ 00), Q{[ct], Pi) ~ [ct](tan/3i — Pi) for i = 1,2, we see that 

J[ct]it) - JfctM = J[ct]M sec A) - Jf,,]([ct] sec/32) 
2 

{cot/3i cos^([ct](tan/3i - - 7r/4) 



7r[ct] 
1 

TTp] 
1 

TTp] 
1 



- cot /?2 cos^([ct](tan P2 ~ P2) - / ^)} 
{cot /5i — cot P2 

+ cot /5i sin(2[ct](tan/3i - pi)) - cot /32 sin(2[ct](tan/52 - P2))} 
[(cot/5i -cot/32){l + sin(2[ct](tan/5i - Pi))} 

+ cot/32{sin(2[ct](tan/5i - Pi)) - sin(2[ct](tan/32 - /32))}] 
-[(cot/?i -cot/32){l + sin(2[ct](tan/3i - pi))} 



cot/52 X 2[ct]{(tan/?i - tan/52) - (/5i - /32)} cos(2[ct](tan/5i - /5i))]. 



Therefore we obtain 



JUit)-JUm 
t-[t] 

r\_/ 

7r(l - c^)f 



{Ci(c)(l + sin(0(c, t))) + t C2{c) cos(0(c, t))} , 



(14) 



12 



where Ci(c) = l/v^T^c^, C2(c) = 2 - (2 - ^/Y^)c^ , and </)(c,t) = 2[ct] (yi^^/c 
— arccosc). For simplicity, we suppose that < a < 6 < 1. Noting Eq. (|T!?|) . we see that 



To get the third inequahty, we used Eq. ((Hj). It is easily obtained that limj^oo -^3(^) = 0, 
since the integral is bounded above by a constant which is independent of t. On the other 
hand, it follows from the Riemann-Lebesgue lemma that lim^^oo ^tSf) = 0. So we have 
limt^oo K2{t) = 0. Thus the proof of Theorem 1 is complete. 

V. CONCLUSION AND DISCUSSION 

In contrast with the classical random walk for which the central limit theorem holds, 
we have shown a weak limit theorem X^'^'' /t dx/ny/l — x"^ {t oo) for the continuous- 
time quantum walk on the line (Theorem 1). Interestingly, although the definition of the 
walk is very different from that of discrete-time one, the limit theorem resembles that 
of discrete-time case, for example, the symmetric Hadamard walk: X^^ /n dx/'n{l — 
x^) Vl — (ri oo). 

Very recently, Romanelli et al. jl^ investigated a continuum time limit for a discrete-time 
quantum walk on Z and obtained the position probability distribution. When the initial 
condition is given by a/(0) = 5/,0)^z(0) = in their notation for the Hadamard walk, the 
distribution at location k and time t becomes the following in our notation: P(^fj^{k,t) = 
J|(t/v^)- More generally, we consider the time evolution given by the following unitary 



[6t]-[at] 



k=0 




dx 



= lim Ksi^t) + lim K^i^t). 



matrix: 




cos 6 sin 6 



sin 6 — cos 6 
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where 9 G (0, 7r/2). Note that 9 = n /A case is equivalent to the Hadamard walk. Then we 
have Pi^{ifi){k,t) = J^{t cos 9). In this case, a similar argument in the proof of Theorem 1 
implies that if — cos6' < a < b < cos 9, then 

P{a<xi^''^^/t<b) [ ^ ^ -.dx, (t^oo), 

J a 'nycos^9~x^ 

where X]. ' denotes a continuous-time quantum walk whose probability distribution is given 
by P(^R0){k,t). As a consequence, we obtain 

-'Vt)"") ^ cos'" ^ X (2m - l)!!/(2m)!!, (t ^ oo). 
In particular, when m = 1, the limit cos'^^/2 is consistent with Eq. (30) in Romanelli et al. 
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